Anomalous quantum reflection of Bose-Einstein condensates from a silicon surface: 

the role of dynamical excitations 
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We investigate the effect of inter-atomic interactions on the quantum-mechanical reflection of 
Bose-Einstein condensates from regions of rapid potential variation. The reflection process depends 
critically on the density and incident velocity of the condensate. For low densities and high velocities, 
the atom cloud has almost the same form before and after reflection. Conversely, at high densities 
and low velocities, the reflection process generates solitons and vortex rings that fragment the con- 
densate. We show that this fragmentation can explain the anomalously low reflection probabilities 
recently measured for low-velocity condensates incident on a silicon surface. 
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Advances in atom cooling and trapping have enabled 
the exploration of diverse quantum transport and inter- 
ference effects for both individual atoms [1] and Bose- 
Einstein condensates (BECs) [2-17]. In BECs, inter- 
atomic interactions can cause complex collective exci- 
tations of the atom cloud, such as soliton and vortex 
formation [8-10], and dynamical and Landau instabil- 
ities [11,12,15-18], which strongly influence the BECs 
center-of-mass motion. There is great current interest in 
using condensed matter systems to provide new tools for 
manipulating ultra-cold atoms, for example micron-scale 
current-carrying wires on a semiconductor surface [19]. 
Understanding and controlling interactions between the 
atoms and the surface is crucial for the development of 
these "atom chips" . Recent experimental studies of Na- 
atom BECs bouncing off a Si surface provide a powerful 
probe of such interactions [13] . When the approach veloc- 
ity, Vx, is low, the rapid variation of the surface potential 
causes quantum-mechanical reflection (QR) of the atom 
cloud. Although quantum-mechanical calculations pre- 
dict that the reflection probability for a single atom will 
increase monotonically as Vx decreases [13,14], the mea- 
sured reflection probability for the BEC decreases with 
decreasing Vx below ~ 1.7 mm s^^: an interesting obser- 
vation that is not yet understood [13]. 

In this Letter, we show that inter- atomic interactions 
can have a pronounced effect on the QR of BECs from Si 
surfaces, consistent with the anomalously low reflection 
probabilities observed in experiment [13]. Our results 
demonstrate that the underlying physics is a generic fea- 
ture of the reflection of BECs from regions of rapid po- 
tential variation, attractive or repulsive. As the BEC 
approaches the surface, its density proflle is modulated 
by the standing wave formed by superposition of incident 
and reflected matter waves. Due to the inter-atomic in- 
teractions in a BEC, this modulation can have a dramatic 
effect on the reflection process. In particular, for high 



density BECs and low incident velocities, it triggers the 
formation of dynamical excitations (solitons and vortex 
rings) that fragment the atom cloud. Our analysis shows 
that although this fragmentation has no intrinsic effect 
on the reflection probability, it disperses the atoms and 
can therefore produce an apparent reduction in the reflec- 
tion probability at low velocities, as observed in experi- 
ment [13]. Our calculations enable us to identify regimes 
in which the fragmentation of the atom cloud that accom- 
panies low-speed QR should be either strongly enhanced 
or suppressed, and thereby propose new experiments to 
test our interpretation of the reflection process. 

We consider BECs containing N ^■^Na atoms in a cylin- 
drically symmetrical harmonic trap. First we take pa- 
rameters from recent experiments [13]: iV = 3 x 10^ 
and longitudinal (radial) trap frequencies Ux (w, ) = 27rx 
3.3 (27r X 4.5) rad s~^, which produce an atom cloud of 
equilibrium peak density no = 2.2 x 10^^ cm~^. At time 
t = ms, we suddenly displace the harmonic trap by a 
distance Ax along the a;-direction [13], and hence accel- 
erate the BEC towards a region of rapid potential energy 
variation at a; = Ax. We consider three different poten- 
tial profiles: potential I [V[ [x), solid curve in Fig. 1(a)] is 
a sharp potential step of height Vg = lO"'^" J; potential 
II [Vjj (x), solid curve in Fig. 1(b)] is an abrupt poten- 
tial drop of depth Vg] potential III [Vjjj (x), solid curve 
in Fig. 1(c)] is a model of the Si surface, whose form is 
specified later. The model potentials I and II enable us 
to identify the key physical processes that occur when a 
BEC undergoes reflection, and thereby develop the un- 
derstanding required to interpret experimental studies of 
BECs reflecting from the more complicated Si surface 
potential. Following the trap displacement, the total po- 
tential energy of each Na atom (mass m) in the BEC is 
Vt (x, r) — V^step (^) + 5"^['^x ~ Ax)^ -f uj'fr'^], where 
r is the radial coordinate and V^gtep {^) is {x), V\\ {x), 
or Viii {x). 
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FIG. 1. Solid curves: schematic representations of Vj (a;) 
(a), Vjj (a:) (b), and Vjjj (x) (c). Dot-dashed curves: the po- 
tential energy of the harmonic trap, shown as a function of 
X, immediately after the trap displacement. Shaded areas in 
(a)-(c) represent the initial atom density profile |$ (a;,0, 0)|^. 
Dashed line in (c) shows the imaginary potential used to 
model the adsorption of Na atoms on the Si surface. 

We determine the dynamics of the BEC in three di- 
mensions by using the Crank-Nicolson [15] method to 
solve the time-dependent Gross-Pitaevskii equation [8,9] 
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where is the Laplacian in cylindrical coordinates, 
a = 2.9 nm is the s-wave scattering length and 4* (x, r, t) 
is the axially symmetrical condensate wavefunction at 
time t, normalized such that l^*! is the number of atoms 
per unit volume. The equilibrium density profile of the 
BEC in the harmonic trap is shown in Fig. 2(a). 

We now explore the dynamics of the BEC in po- 
tential I. When the trap is displaced at t = ms, 
the initial potential energy of the BEC is increased by 
AV ~ ^moj'^Ax'^ , as shown in Fig. 1(a), causing the 
atoms to be incident on the potential step with a veloc- 
ity Vx ~ LOxAx. For the displacements considered here, 
AV << Vg. Hence, in a classical picture, we expect the 
BEC to be completely reflected by the step. Figure 2 
shows images of the BEC undergoing reflection from po- 
tential I after a trap displacement of 60 fim, for which 
Vx ~ 1.2 mm s"^. At i = 90 ms [Fig. 2(b)], approx- 
imately half the atoms have reached the barrier, and a 
standing wave has formed as a result of the superposition 
of the incident and reflected waves, causing periodic mod- 
ulations in the density profile. Due to the inter-atomic 
interactions, the high density at the peaks of the stand- 
ing wave causes atoms to be pushed out from the axis of 
cylindrical symmetry (the x-axis), thus transferring mo- 
mentum into the transverse direction. These atoms ap- 



pear as the arrowed jet-like "lobes" at the top and bottom 
of Fig. 2(b). The lobes ride up the walls of the harmonic 
trap and are then reflected back towards the x-axis. This 
causes a cylindrically symmetrical soliton to form, which 
appears as the two white stripes (arrowed) in Fig. 2(c). 
The soliton then decays via the snake instability [20] into 
two vortex rings centered on r = 0, which appear as pairs 
of density nodes within the dashed box in Fig. 2(d). The 
phase of the BEC wavefunction within the dashed box 
in Fig. 2(d) is shown to the right of the density profile. 
The figure reveals a 27r phase change around each den- 
sity node, indicating quantized circulation. At the end of 
the oscillation, the atom cloud still contains vortex rings, 
and has a fragmented appearance [Fig. 3(a)]. Similar 
dynamics have been observed for two-dimensional BECs 
reflecting from a circular hard wall [10]. 
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FIG. 2. Grey-scale plots of atom density (black — high) in 
the X — r plane (axes inset) for a BEC with no — 2.2 x 10^^ 
cm"^, ujx (t^r) = 27r X 3.3 (27r x 4.5) rad s"""^ and ~ 1.2 
mm s~^, reflecting from potential I, shown at t — ms (a), 
90 ms (b), 122 ms (c) and 143 ms (d). The phase of the BEC 
wavefunction within the dashed box in (d) is shown in the 
grey-scale plot (white = 0, black — 2n) to the right of the 
density profile. Horizontal bar shows scale. 

For disruption to occur, the time taken for lobes to 
form, w Ir/^s, where Ir is the BECs radial width 
and Vs ~ ^jK^nQa/ir-m? is the average speed of sound 
[2], must be less than the reflection time Ir k, Ix/vx, 
where Ix is the BECs longitudinal width. We therefore 
estimate that reflection will not significantly disrupt the 
BEC when 

Vx ^ Vs (Ix/lr) , (2) 

or, equivalently, when Aa; {vslx) / {Ir'-^x)- Our numeri- 
cal simulations confirm this [21]: for a large Vx ~ 2.1 mm 
s"^ {Ax — 100 /im), the reflected atom cloud aX t — 150 
ms [Fig. 3(b)] contains no excitations and has a smooth 
density profile, similar to the original BEC ground state 
[Fig. 2(a)]. Similar behavior has been reported for BECs 
on Bragg reflection in optical lattices [9,16,17,22]. 

We now consider QR of a BEC from the abrupt poten- 
tial drop (potential II) shown in Fig. 1 (b) . The reflected 
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atom clouds at t = 150 ms for « 1.2 and 2.1 mm 
s~-^ are shown in Figs. 3(c) and (d) respectively. They 
are qualitatively remarkably similar to the equivalent im- 
ages for potential I [Figs. 3(a) and (b)], except that they 
are slightly smaller because fewer atoms are reflected. 
When Vx ~ 1.2 mm s~^, the reflected atom cloud has a 
fragmented appearance and contains a vortex ring, en- 
closed by the arrows in Fig. 3(c). By contrast, when 
Vx ~ 2.1 mm s^^ the BEC contains no excitations and 
has a smooth density profile [Fig. 3(d)]. This is because, 
just as for potential I, when Vx ~ 2.1 mm the reflec- 
tion process is too fast for side lobes to form in response 
to the standing wave. 
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FIG. 3. Grey-scale plots of atom density (black — high) in 
the x—r plane (axes inset) for a BEC with with no = 2.2 x 10^^ 
cm~^, ujx {<-0r) — 2-kx 3.3 (2n x 4.5) rad s~^ and Vx ~ 1.2 mm 
s~^ (a), 2.1 mm (b), reflecting from potential I, shown at 
t = 150 ms. (c) & (d) Same as (a) & (b), but for potential II. 
(e) & (f) Same as (a) & (b), but for potential III. Horizon- 
tal bar shows scale. Arrows in (a), (c) and (e) indicate the 
direction of quantized circulation around vortex cores. 

Having studied two simple abrupt potential steps, we 
now investigate QR from a Si surface [13] at x = Ax, 
which we model by potential III [Vjjj {x), solid curve 
in Fig. 1(c)]. Potential III is based on the Casimir- 
Polder potential Vcp (x') = -Ci/x'^ {x' + 3Xa/2Tr'^) , 
where d 9.1 x 10"^^ J m"*, Aa = 590 nm is the effec- 
tive atomic transition wavelength [23] , and x' = x ~ Ax 
is the distance from the surface. For x < Ax — S, 
where 6 = 0.15 fim is a small offset from the surface 
[Fig. 1(c)], we set Vjjj (x) — Vcp {x'). However, since 
Vcp — > — CX3 as a; — > Ax, for x > Ax — S we set Vjjj (x) = 
Vcp {S) + i{x - Ax + S) V^^, where V^^ = 1.6 x lO'^^ J 
m~^. No atoms are reflected when x > Ax — 5 because 
the real part of (x) is constant. The imaginary part 
of Viii (x) [dashed curve in Fig. 1(c)] models the adsorp- 
tion of Na atoms by the Si surface [24] . We expect that 
the BEC will undergo QR when \dX/dx\ > 1 [14], where 
A is the local de Broglie wavelength. For Vx = I mm 
s~^, the atoms are most likely to reflect ~ 1/ini from the 
surface, where \dX/dx\ is maximal [25]. 

Figures 3(e) and (f) show the reflected atom clouds at 



t = 150 ms for Vx ~ 1.2 and 2.1 mm s""'^ respectively. Al- 
though they are smaller than the corresponding images 
for potentials I and II because fewer atoms are reflected, 
the qualitative behavior of the BEC is the same. When 
Vx ~ 1.2 mm s~^, the atom cloud is disrupted and con- 
tains a vortex ring [Fig. 3(e)], whereas for Vx ~ 2.1 mm 
s~^ the BEC has a smooth density profile and contains no 
excitations [Fig. 3(f)]. This is consistent with the exper- 
iments of Pasquini et al. [13], who observed "excited and 
sometimes fragmented" atom clouds following the QR of 
BECs incident on a Si surface at low Vx- 

To quantify how the QR process depends on Vx, and 
to relate our simulations to experiment [13], we made 
numerical calculations of the probability, R, that each 
atom in the BEC will quantum reflect from the Si sur- 
face potential. The solid curves in Fig. 4 show R versus 
Vx calculated for atom clouds with (squares) and with- 
out (crosses) inter-atomic interactions. In each case, R 
increases with decreasing Vx- Similar behavior was re- 
ported for single atoms in Refs. [13,14,26]. 
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FIG. 4. Solid curves: R{vx) for potential III, calculated 
from the entire reflected atom cloud with (squares) and with- 
out (crosses) inter-atomic interactions. Dotted curve with 
circles: Riv^) calculated with inter-atomic interactions, but 
omitting regions of the atom cloud where |^[^ < 0.25no. 
For the solid and dotted curves, no — 2.2 x 10^^ cm~^ 
and LOx (<^r) = 27r x 3.3 {2-k x 4.5) rad s~^. Dot-dashed 
curve with triangles: as dotted curve with circles, but for 
a BEC with no = 3.6 x 10^^ cm~'^. Inset: data points 
(squares) show In («p) vs. In (no), where Vp is the position of 
the peak in R (vx ) curves calculated omitting regions where 
1^1^ < 0.25no, and dashed line is least squares fit. 

A key feature of Fig. 4 is that the R{vx) curves cal- 
culated with and without inter-atomic interactions are 
nearly identical [27], indicating that the interactions have 
little intrinsic effect on the reflection probability [28]. 
Consequently, inter-atomic interactions do not immedi- 
ately explain why experimentally measured R values de- 
crease with decreasing Vx ^ 1.7 mm [13]. However, 
a possible explanation of this anomalous behavior does 
emerge when we consider how the interactions affect the 
structure of the BEC at low Vx- In the absence of in- 
teractions, QR does not change the form of the atom 
cloud. By contrast, interactions cause disruption and 
fragmentation of the atom cloud when Vx ^ 1.7 mm s^^ 
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{Ax < 80 /Ltm): precisely the regime where anomalously 
low R values are observed in experiment [13]. Severe dis- 
ruption of the cloud can make detection of atoms difficult 
in low density regions of the BEC [29] , facilitate thermal 
excitation, or even eject atoms from the trap [12], thus 
accounting for the anomalously low R values observed in 
experiment. To simulate the effect of such losses, we re- 
calculated R omitting the contribution from atoms in low 
density regions of the BEC for which l^*]^ < ric, where Uc 
is a threshold below which detection of atoms can be dif- 
ficult [29]. The dotted curve with circles in Fig. 4 shows 
R {vx) calculated for Uc = 0.25 uq. In this case, as in ex- 
periment [13], R decreases with decreasing < Vp ~ 1.7 
mm s~"'^; in the regime where fragmentation of the atom 
cloud produces large areas of low density [30] . 




disrupted by QR [Fig. 3(f)]. 

In summary, we have investigated the reflection of 
BECs from three abrupt potential steps, and shown that 
the dynamics are qualitatively the same in each case. 
For high Vx, reflection occurs cleanly, causing almost no 
change in the atom density profile. But at low v^, the 
standing wave formed from the superposition of the inci- 
dent and reflected waves generates dynamical excitations 
that disrupt the atom cloud. This disruption could ex- 
plain the anomalously low R values recently measured 
for BECs incident on a Si surface at low velocity. Frag- 
mentation of the atom cloud is particularly dramatic for 
dense slow-moving cigar-shaped BECs, and further ex- 
periments in this regime would test our interpretation of 
the reflection process. Conversely, for low-density BECs, 
which are not significantly disrupted by QR, the observed 
values of R should continue to increase with decreasing 
Vxi even when is small. In this regime, QR near a Si 
surface should provide an effective means of containing 
BECs and decoupling them from the surface. 

We thank T. Pasquini and O. Morsch for helpful dis- 
cussions. This work is funded by EPSRC UK and ARC. 
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FIG. 5. Grey-scale plots of atom density (black = high) in 
the X — r plane (axes inset) for a BEC with no = 2.2 x 10^^ 
cm~^, ujx {ijJr) = 27r {2-k x 4.5) rad s"""^ and Vx ~ 2.1 mm , 
reflecting from potential III, shown at f = ms (a), 270 ms 

(b) , and 460 ms (c). Horizontal bar shows scale. Arrows in 

(c) enclose a vortex ring. 

To explore the parameter space, we repeated our cal- 
culations for BECs in the same harmonic trap as before, 
but with different values of N and tiq. When we increase 
N to 10^ (no = 3.6 x 10^^ cm'^), the peak in our cal- 
culated R{vx) curve (dot-dashed in Fig. 4) occurs at a 
higher Vp of 1.9 mm s~^ (corresponding to Ax = 90 /im). 
Simulations for BECs with a range of N reveal that a plot 
of In (vp) versus In (uq) (Fig. 4 inset) is linear, with a gra- 
dient of 0.5 ± 0.01. This is because the peak in R{vx) 
coincides with the onset of disruption of the atom cloud 
so that, from inequality (2), Vp cx y/no (Ix/lr) oc y^rio: 
(Ix/lr) is independent of N and no for the trap considered 
so far. According to this analysis, for given no, Vp will be 
higher for a more cigar-shaped BEC. To confirm this, we 
investigated QR for an elongated atom cloud [Fig. 5(a)] 
formed by reducing lUx to 27r rad and increasing N to 
9.8 X 10^ to maintain the original value of no = 2.2 x 10^^ 
cm^'^. For Vx ~ 2.1 mm (Ax = 330 /xm), large side 
lobes form during QR [Fig. 5(b)], leading to fragmenta- 
tion of the atom cloud and production of a vortex ring 
[Fig. 5(c)] [31]. By contrast, for the same values of no 
and Vx, the more spherical BEC shown in Fig. 2(a) is not 
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